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Abstract. In this work we generalise previous results connecting (rational) Gaudin 
magnet models and classical separation of variables. It is shown that the connection 
persists for the case of linear r-matrix algebra which corresponds to the trigonometric 
4x4 r-matrix (of the XXZ type). We comment also on the corresponding problem 
for the elliptic (XYZ) r-matrix. A prescription for obtaining integrable systems 
associated with multiple poles of an L-operator is given. 



1. Introduction. Separation of variables for the Hamilton- Jacobi and Schrodinger 
equations have long been known as methods for explicit solution of these equations 
in appropriate circumstances. The technical requirements for this method of so- 
lution have quite fully developed in recent years (see [1-8]). In particular the 
relationship between the separable systems and the Gaudin magnet [4,9] integrable 
systems models has been established via r-matrix algebra methods, where the r- 
matrix corresponds to the rational or so called XXX case, [4-8]. This relationship 
works very clearly with separable coordinate systems on spaces of constant curva- 
ture. The question we answer here is how these notions can be extended to include 
the so-called trigonometric r-matrix algebra in the XXZ case. To do this let us 
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recall the fundamental ideas of the r-matrix formalism (see [10,11] and references 
in there). For a classical mechanical system the basic (linear) r-matrix algebra is 

(1.1) {L{u) (g) /, / (g) L{v)} = [r{u - v), L{u) ® I ^- 1 ® L{v)] , 

where {-, •} is the Poisson bracket and [-, ■] the matrix commutator bracket. The 
operator L{u) is taken to be the 2x2 matrix 

L{u)ii = -L(w)22 = A{u), L{u)i2 = -B(w), L{u)2i = C(w), 

and r{u) is a suitable 4x4 matrix of scalars solving the classical Yang-Baxter 
equation [10,11]; u being arbitrary constant is called the spectral parameter. In the 
case of the XXZ r-matrix algebra the non zero elements of r can be taken to be 

(1.2) r(«)ii = r(u)44 = coth(«), r(w)23 = r(w)32 



sinh(tt) 

In component form, the r-matrix algebra relations are 

{A{u), A{v)} = {B{u), B{v)} = {C{u), Civ)} = , 
(1.3) {A{u), B{v)} = -^T-^^ — (cosh(u - v)B{v) - B{u)) , 

{-cosh{u-v)C{v) + C{u)), 



1 


sinh(tt — 


v) 


1 




sinh(w — 


v) 


-2 




sinh(« — 


v) 



If we now make the ansatz A{u) = coth(tt)/S'3, B{u) = (1/ sinh(tt))S'+ and 
C{u) = (1/ sinh(w))<S'_ these relations imply 

(1.4) {^3,^±} = ±S'±, {S+,S-} = 2Ss. 

To relate this observation to the separation of variables methods, we form the 
L{u) operator with elements 

^ ^ sinh(i(, — Crv) ^ ^ sinh('U — €„) " ' 

a=l a=l 

(1.5) 

n 

A{u) = ^ COth(w - Cq) Ssa , 

where 

(1-6) {Ssa, S±i3} ^ ±5af3S±a , {S+a, S-p} ^ 25ai3S3a ■ 

The r-matrix algebra relations, (1.1) or (1.3), imply 
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i.e., that det L{u) is a generating function of the constants of the motion. In 
particular we have 

Cry 



det L{u) = A'^iu) + B{u)C{u) = ^ 



a=l 



sinh^(tt — Ca) 



+ Ha COth('U — Cq) 



where Ca = (Ssa)'^ + S+aS-ce are the Casimir elements of the algebra generated 
by elements S+a, S-a and S^a- Furthermore Hq = ^^<^ 

(1.9) Ha=^ (^SsaSsfS COth(ea - 6)3) + . , _ r{S+aS-i3 + S+pS-a)) ■ 

With the following realization of the algebra in terms of the canonical coordinates 

Xa and Pa, {pa,Xf3} = -5q,/3 : 



(1.10) 

the constants (1.9) have the form 



(1.12) Ha^W 



^ t^a sinh(ea - e^) 



{xl^Pp + X^pI, - 2XaXi3PaPi3 COsh(ea - 6/3)) , 



and Hq = XaPa ■ Notice that all Cq = in such a representation. 

2. Variable Separation for the XXZ Magnet. Proceeding as in [4,7,8], we 

choose separable coordinates such that B{u) = 0, i.e., u = uj, j = 1, . . . , n — 1. 
This implies 5^^x^/sinh(M — Cq.) = for tt = ui, . . . , Un-i, which in turn implies 
that we choose coordinates according to 



n^=i'sinh(«,-ea) 



(2.1) 



" n^5^aSinh(e/3 - Co,) ' 

motivated by the general formula 



(2.2) 



sinh(« — Cq;) n^=i sinh(tt — e^) 



For each Uj we can define the canonically conjugate coordinate Vj as follows: 

1 

(2.3) Vj = A{uj) = coth{uj - e^) XaPa , ^ < j <n - 1 , Vn = Hq . 



Oi=l 



The coordinates ttj, Vj {i,j = l,...,n) satisfy the canonical bracket relations 
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The changing of variables Xa,Pa for the new variables Wj, Vi is the procedure of 
separation of variables. The matrix elements of the L-operator can be expressed in 
terms of these variables according to the formulas 



A{u) =iB{u) [2 cosh [u + TT^zluj - T^l^^e^) 



(2.5) 



n-1 



-2vj njj^i sinh(w^ 



sinh(u — Uj) II/j^j sinh(uj — Uk) 
The entry C{u) can be computed by using the formula 



Pa =Xae-''- [2cosh(Sp^ 



Un 



—2vj ^a=i sinh(wj — Ca) 
^ sinh(ea - Uj) U^^j sinh(ttj - Uk) 



n-1 



(2.6) 



+ E 



This gives the relation between the coordinates Xa,Pa and Ui, Vi where = Hq. 
The equation for the eigenvalue curve F : det(L(u) — XI) = 0, has the form 

- A{u)^ - B{u)C{u) = . 

If we put u — Uj, j — 1, . . . ,n — l into this equation then A = ±fj. Thus variables 
Uj and Vj (j = 1, . . . ,n — 1) lie on the curve T: 



(2.7) 



Vj — Ha coth(ttj — Cq,) — Hq =Vj+ det L{uj) = , 



a=l 



Equations (2.7) are the separation equations for the degrees of freedom connected 
with the values of the integrals H^- (Note that Yla=i = 0.) 

For illustrative reasons it is more transparent to use the variables Ai = e^"- and 
Ui = e"*. Then many of the expressions given have algebraic form. For example, 
the nearest object we have to a Hamiltonian in the case of XXZ r-matrix algebra 
is H = Yl7=i ^l^i which has the form 



n-1 



i=l 



n-1 



(2.8) = {4U]^,A,)-\-ir (pI. + Hi) ^ 

Note that 



Ilj^iSmh.[ui — Uj) 

ni^Ml-uf) 



i=l 



_ Ui--- Un-1 sr^ A 2 

e — , , ^aXa ■ 



Ai---An 



a=l 



If we adopt the standard procedure and write pui — dW/dUi then the equation 



H = E admits separation of variables via the usual ansatz W = Yl^=i Wi{Ui). The 
separation equations can be written in the alternate form 



[uUiuf-AD] 



BTT. 



n-1 



2fc-2 
3 ' 
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where the A/, are related to the Hj via 



ai<...<a„-fc ^ ' 



An-1 = E — Hq ^ . 



fc=l 

From what has been developed so far we see that separation of variables goes 
through for XXZ r-matrix algebras constructed in this way. In the previous article 
[8] for the case of spaces of constant curvature we essentially have the rational r- 
matrix algebra and it is possible to formulate using well defined limiting procedures 
the cases of integrable systems for which some of the parameters are equal. What 
was also established previously was the construction of integrable systems given on 
the algebra with commutation relations 

(2-9) {{Zj)e, {Zi )m} = -Sjj'{Zj_^f,_j^j) seems, 

where £,m,s= 1, 2, 3, < j < Nj, < k < Nj', < J < p and e^ms is the usual 
totally antisymmetric tensor, and the vector has the form 

(i 12iiPiPj + l-i + ^i^j + l-i) 
2 Z^iPi -^j+l-i 

in the coordinate representation. Indeed, if we adopt the limiting procedure 
A-^ ^ A( + -^e]_^, j = l,...,Nj, J=l,...,p, 

Nj 

E 

i=2 

Nj 



\ i=2 

/ Nj 

(2.11) v«/hi+E'^5+i-i^ 



where 



i=2 



JJ-1 _ jii (J A _JA\ fjJ - 1 



and A'^i + ... + Np = n + 1, then the Hamiltonian H has the form 
(2.12) H = (n»«^,)-' 5: + Hi) 

with obvious separation equations. (We require "^eg = and take the limit as the 

J,l . n C U 1 AT 1 rol ^ rpi „ „ . • , r, • t „ J- J-„ 
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can be derived by applying these procedures to det L('u). We will, however, adopt 
a different and more general strategy. If we leave the matrix elements of L{u) in 
the form (1.5) and subject the resulting expression for det L{U) 



ldetL{U) = ^(^j2Ss^^ + 



2U' 



EE 



i=l \a=l^^ ^« 



(2.13) 



S: 



3a 



a—l " 



S: 



3a 



1 U^-Al 



where U = e'^, and S±a = Sia ± iS2ai to the transformations 

M + i = i,...,iVj, j = i,...,p, 

J 



Nj 



then we arrive at a general expression for the generating function det L(C7). The 
constants of the motion are obtained by the usual means of expanding the ex- 
pression following from (2.13) in partial fractions and reading off the independant 
components. In the case of degenerate roots the expression can be readily modified. 
Accordingly we have 



\dei L{U) = \[j^{Zi,^))l 

E EE^ 



+ 2C/^ 



i=\ \J=1 i=0 
p Nj-\ ^ I ^ a \ J 



9 V A, \ J 



(2.14) 



d 



dAj U^-A] 



From this expression constants of the motion can be deduced just as before. The 
separation of variables proceeds as usual in the case of the choice of coordinates as 
given in [8]. The expressions for the coordinates corresponding to multiple roots 
with signature Ni, N2, Np can be obtained from the generic case by the limiting 
procedures already outlined. In rational form the generic coordinates are 

(2.15) = [U]^,{A'j - U!m:=lUJ ^'tf\2s ■ 

For the case of signature Ni, N2., Np the coordinates are given by the relations 

'J -2 



lr=l 

/^AT,\ A T * \ f-\ A T I A T-n A T.\ AT, I * 



i=l 
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This gives a complete description of the separation of variables procedure for the 
signature A^i, A^2, -/Vp case. We illustrate these ideas with two examples. 

A. The case of signature 2,1 and dimension 3. In this case the generating 
function assumes the form 

det L{u) = —-^ iUzlZl + ZlZl) - / mUZl), 

smh [u — ei) smn^ei — 63) 

+(Zi)i(Zf)i) +cosh(ei -e3)(Z!)3(^i')3} 

2smh (u - ei) L smh(ei - 63) ' 

+ cosh(ei - e^){ZlUZl)^ - —-j^ {{ZlUZl)^ 

smh (ei — 63) 

+ cosh(ei - eMZl)2{Zl)2 + {Zl),{Zl),))] 

- coth(« - 63) \ ^ (- {{zMzD, + {zMzl),) 

y smni^ei — e^,) 

+ cosh(ei - e3)(Zf)3(Zi)3) + —-jj^ {{ZlUZl)^ 

smh (ei — 63) 

+ cosh(ei - 63) {{Zl)^{Zl)^ + {Zl)2{Zl)2))]. 

The constants of the motion can be deduced from the coefficents of independent 
functions of u. In the coordinate representation these constants have the form 

Hi = x\pI + p\xl - 2X1X2P1P2 + . ,1 ^ r + Avi - 2a;i3;3PiP3) , 

smh(ei — 63) ^ ' 

2 

Hi = ^—rj T ixiX2pl + PiP2xl - {X1X3P2P3 + P1P3X2X3) cosh(ei - 63)) 

smh(ei — 63) ^ ^ 

(xips+pixsf, 



sinh (ei — 63) 



where we have used the notation xi = x\ , X2 = X2 and x^ = xf , with similar 
relations for the pi's. The coordinates are given by the formulas 



X-, = 



sinh(wi — ei) sinh('U2 ~ ^1) 



sinh(ei — 63) ' 
sinh('Ui — ei) sinh(w2 — ei) cosh(ei — 63) 1 

sinh^(ei - 63) sinh(ei - 63) 

X (sinh(iii — ei) cosh(t(,2 — ei) + sinh(ui — 63) cosh(u2 — 63)) , 
sinh('Ui — 63) sinh('U2 — 63) 



^3 



sinh (ei — 63) 
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has the form 

+ (^2,^2)5 16-^1 [5^1-^1 + (^1)2(^2)2 + 2A, ((^1)3(^2)3 + (^1)1(^2)1)] 



+ (^723^2^4 [33Zi.Zi + 36Z1.Z2 + SAjZi.Zs + ^A\Z^.Z^\ 

+ ^^7^4:42)3 2^? [26A1Z1.Z2 + %A\Z^.Z^ 
+UAIZ1.Z3 + Zi.Zi + 4AIZ3Z2 + {Zi)l] 
+ (^2,^2)2 2 [AiZs.Zs + 6AIZ2.Z, 

+ Aj (5 ((Z2)i + (Z2)?) + 6(^3)2(^1)1 + 4(^2)2) 

+Ai (6(Z2)2(Zi)2 + (Z2)l(Zi)i + 4(Z2)3(Zi)3 + (Zif)] 

+ ^2^^2 2 [^?^3.^3 + 2^iZ2.Z3 + iZ2)l + (^2)? + (^1)3(^3)3] • 

The coordinates are given by 



4 ' ^ 4 Vt^2 C/i/ 4t/it/2 ' 



8^1 8 \U2 UiJ 8U1U2 ' 



8^f 8 Vt/2 C^i/ 8t/it/2 ■ 

3. The XYZ Magnet. These methods can be extended to the case of eUiptic or 
XYZ r-matrix algebras. The only difference is that in this case a solution of the 
problem via separation of variables is not yet known^ but the coalescing of indices 
goes through just as before. Indeed, the operator L{u) can be taken just as in (1.1). 
The non zero elements of the r-matrix in this case are 

cn('u) 1 - dn(w) 

r(tt)ii = r(u)44 = — — , r('u)i4 = r('u)4i - 



(3.1) r(tt)23 = r-(tt)32 



sn(w) 2sn(w) 
1 + dn('u) 



2sn(w) 

We now make the ansatz 

A{u) = ^^3 , B{u) = [(1 + dn(«))5_ + (1 - dn{u))S+] , 

(3.2) 

^^""^ ^ 2^) + Mu))S+ + (1 - dn(«))5_] . 



■"■See [12] where the variable separation has been done for the periodic classical XYZ-chain 
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Here the S±, S3 obey the same commutation relations as (1.4). We choose the L{u) 
operator to be 

(3.3) 

^ cn{u-ea) 
f-^ sn(« - ea) 

n ^ 

B{u) = ^ - — ^— r [(1 + dn{u - ea))S'_« + (1 - dn(w - ea))5'+«] , 

a=l ^^^l'" ^a/ 

n 

'^(") = E o / _ ^ [(1 + dn(w - e«))^+a + (1 - dn{u - e,))5_,] . 

a=l ^^^\^ ^a/ 

The determinant of L{u) is once again a generator of the constants of the motion. 
It has the form 

n 

(3.4) det L{u) = HaE{u - + iK') + Hq 

a=l 

where 

Ha = ^ ^- [SiaSii3 + dn(ea - ei3)S2aS2p + Cn(ea - e^)S'3a'S'3/3] , 

i^o = 2A;2^ ' E^^l^-^ [S^aSip + dn(e« - ep)S2aS2p + cn(e« - e^)53a^3/3] 

a, /J ^ ' 

(3.5) 



- ' [k'^cn{ea - ep)S2aS2p + dn(e« - ep)S3aSzp\ - Y^ik'^Sl^ + S[ 

a, 13 a=l 



3a] 



Here E[z) = dn^{u)du is Jacobi's epsilon function. The same is now true as 
for the case of XXZ r-matrix algebras: if we subject the e^'s and the S'j^'s to the 
transformations given by (2.11), then we arrive at the generating function for the 
constants of motion for a root structure having the signature Ni,N2, ....,Np. The 
expression for this function is 

3 /iVj-1 / a \r \^ 

(3.6) detL{u) = Yij2 [dJ^J " 1 ' 

where fi{z) = l/sn(2), f2{z) = dn{z)/sn{z), and /a (2) = cn(2;)/sn(2;). 

As an example, the generating function corresponding to signature 2,1 is 

det L{u) = HiE{u - ei + iK') + H2E{u -63 + iK') + Hs + —rj^ rH^ 

sn*(tt — d) 

cn(^-el)dn(^-el) \ 1 

H 57 ^ -"5 H ^7 T-ne H ^7 ^-"7, 
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where 



/c^cn(ei - e3)dn(ei - 63) ^ . , cMei-ea) 2 1 
sn(ei - 63) sn(ei - 63) 

- 2fc2sn(ei - e3)(Z2)i(Z|)i + 2dn(ei - e^){Zl)^{Z\)^ , 
E./^ — Z\.Z\, }1^^ = 2Z\.Z\^ 

He = z\.z\ - {zDl - {zDl - k\zl)l + A zluzl), 

sm^ei 63^ 



sn(ei - 63) 



We note that the ideas developed here also work in the case of separation of 
variables for spaces of constant Riemannian curvature, as developed in previous 
articles [6-8]. Indeed, in that case the rational r-matrix algebra is as before and 
the non zero elements of the r-matrix are 

(3.7) r(tt)ii = r{u)44 = r{u)23 = r{u)32 = 1- 

The generating function of the constants of the motion for signature A^i, ....,Np is 
then 

3 / P Nj-l fyj \ 




(3.8) detL(«)=j: j:^:^-^ 

k=l \J=1 j=0 ^ ^' 

This is the generalisation of the generating function for separable coordinates on 
spaces of constant curvature of dimension n = Y7j=i + 1- Indeed, if we use the 
form (3.8) and if e/. = for k = 1,2,3 then we have the generating function on 
the sphere for generic ellipsoidal coordinates, and if ei = — 1/4, 62 = 1/4, 63 = 
then we have the generating function of ellipsoidal coordinates in n-dimensional 
Euclidean space. As an example consider the system with signature 2,1. The 
generating function is then 

(£fM£i) _ (ziM£i) (zluzl) _ (zfuzl) 

det L{u) = i^izf^ + 

u — ei [u — eiY 

{u — ei)^ {u — ei)^ w — 63 {u — 63)^ 
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